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ABSTRACT 

i 

Cases simulating beam currents to about 29 MA were done. Significant noise reduc- 
tion and code speed-up were achieved for these cases. The net current scaling observed at 
‘ lower beam current ( I nt t ~ Iblam) continued to hold up well, as did the equilibrium itself. 

A slight deviation in the expected thickness and radial structure of the current layer was 
observed, however. A laminar flow model is developed that appears to give good agreement 
with the simulation results. Suggestions for future work are discussed briefly. 
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1.0 INTRODUCTION 



Using our earlier test data (reproduced in Appendix A), in order to reduce noise in 
the high-current cases we increased the number of particles per cell to 375 for the sheath, 
used radial current smoothing over three adjacent cells and eliminated current inside the 
axial cell ( R ~ 0) to reduce the fictitious spike near the axis. This decreased the fluctuation 
amplitude by about a factor of three and made estimation of the average current values 
from the plots much easier. 

The summary plots shown here are in the same format as those in Appendix A but 
continue the beam-current analyzed on to much higher levels. The trends are essentially 
as expected. 

We have completed a laminar-flow equilibrium model that predicts the radial profile 
of the axial velocity of electrons and the relationship between azimuthal magnetic field 
and radial electric field. Specific predictions are obtained for the net current and ion-beam 
temperature necessary for equilibrium. In all cases the agreement between the model 
predictions and the ISIS simulations is excellent. The model also predicts parameter 
scalings which have not yet been tested in the simulations and, therefore, need further 
verification. 




2.0 RESULTS OF 1.8-28.8 MA SIMULATIONS 



Table 1 describes the simulation parameters and selected results for the new runs 
carried out in August 1989. We have also Included a few older simulations that are used 
in some of the summary plots. In all these runs, the axial velocity of the ion beam was 
Pai = ,0.62, the annulus between the beam and the conducting pipe wall was filled with 
plasma at five times the beam density, and the ion beam temperature was T, = 43.2 kV 
(the same as concluded at the end of the research in Ref. 1). 

Cases were run at 1.8 MA for both a normal configuration and one with the uninter- 
esting inner region of the beam, which usually shows (nearly) zero net current, replaced 
with a particle reflecting metal boundary at R = 25 cm (see Figs. 1-5). These cases 
agreed within ~ 10 percent, so higher currents were run with this feature which decreases 
the simulations time by about 40 percent. Because the CPU run time to steady state for 
the net current increases roughly as the square root of the beam current, we were able 
to run the 7.2 MA and 28.8 MA cases readily (Figs. 6-9 and 10-13, respectively) but a 
115 MA case (which might require 15-20 CPU hours) was not undertaken due to the cost. 
The width of the nonzero net current region, which had previously scaled above 100 kA 
as approximately ~ continues its rapid decrease for h tam > 7 MA, (Figs. 4, 5, 9, 

and 13). The radial current profiles of these and some lower current cases run previously 
are superimposed in Fig. 14. All cases continued to exhibit approximately the same net 
current scaling (I net ~ 15-17) as observed at lower currents. 

The radial structure of the sheath region is more complex than earlier thought. Noise 
in earlier runs masked this structure. A negative current region at the inside edge of 
the sheath (e.g., Fig. 14, 1.8-18.8 MA) became more obvious at higher currents. This 
region is a local over-neutralization of the current at the edge of the beam electrons. The 
sheath electron cyclotron orbit about Bt imparts a large axial velocity component. At its 
innermost radial position there is some overlap with original beam electrons which have 
not yet been lost to the wall (see Fig. 21). AI 30 , the sheath electrons at their inner radial 
position overneutralize the beam charge, form a virtual cathode (reverses the sign of E r ), 
and loose all their energy (Figs. 21a, c, and e). In this “mixing” region the sheath electrons 
lose energy and axial momentum and the original beam electrons gain it (see Fig. 21d-f). 




The complicated phase space of the beam electrons (Figs. 21b, d, and f) as they leave 
radially is due to the field reversal of E r and B$ in the sheath region just Inside the beam 
edge. E r and Bf have very similar radial profiles and go to zero at the same radial position 
(Figs. 1 and 2, 6 and 7, and 10 and 11). 

No significant compression or expansion of the ion beam was observed in any of these 
runs, so we postulate that the equilibrium ion-beam temperature is relatively independent 
of current. 

Figures 18-20 show that a steady state has been reached at least in the net current. 
The noise has been substantially reduced over the level observed in the earlier simulations 
reported on in Appendix A. 
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TABLE 1. SIMULATION RESULTS AND PARAMETERS. 



LLNL Plasmoid Propagation Results August 1989 







Cell 






Number 


Run 




Size 
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1.8 MA 


0.1 


-0.05 


53.25 
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0.1 


25 


51.8 




AC 


7.2 MA 
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25 


50.25 




AD 


28.8 MA 


0.025 


25 


50.25 




AL 


22.5 kA 


Run done at AFWL Aug-Oct 1988 


BE 


452 kA 


Run done at AFWL Aug-Oct 1988 


BB 


1.8 mA 


Run done at AFWL Aug-Oct 1988 
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(b) 



Figure 1. Periodic simulation results comparing azimuthal magnetic field B$ (B3, 0.6 kG) 
versus radial position X2 = R (X2, cm) at ct = 400 cm for two runs with 
I b = 1.8 MA. (a) run AA, (b) run AB (hollow beam), NOTE CHANGE IN X2 
SCALING. 
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Figure 2. Periodic simulation results comparing E r (E2, 0.51 MeV/cm) versus R (X2, cm) 
at ct=400 cm for two runs with h = 1.8 MA. (a) run AA, (b) run AB (hollow 
beam), NOTE CHANGE IN X2 SCALING. 
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Figure 3. Periodic simulation results comparing axially-integrated J x versus R (X2, cm) 
at ct = 400 cm for two runs with Jj = 1.8 MA. (a) run A A, (b) run AB (hollow 
beam), NOTE CHANGE IN X2 SCALING. 
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Figure 5. Time averaged c 
I h = 1.8 MA. 
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Figure 7. Periodic simulation AC, E, (E2, 0.51 MeV/cm) versus R (X2, cm) at cl 
500 cm for /(, = 7.2 MA« 




J1-IMT1 





LLNL PLASMOID PROPAGATION 



o 




(Yl) INdddOD 



13 



Figure 9. Time averaged current profile for run AC near end of run with I $ 7.2 MA. 
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Figure 11. Periodic simulation AD, E r (E2, 0.51 MeV/cm) versus R (X2, cm) at ct 
200 cm for L - 28.8 MA. 
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Figure 12. Periodic simulation AD, axially-integrated J t versus r at ct = 200 cm for 
I h = 28.8 MA. 
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Figure 13. Time averaged current profile for run AD near end of run with I b — 
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TIME- 400.20000 



jmparing time history and frequency spectrum 
ns. /» = 1.8 MA. Current in units of 1.35 kA, 
x 10 10 rad/s. (a) run AA, (b) run AB (hollow 
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Figure 19. Periodic simulation AC, time history and frequency spectrum of the net cur- 
rent. h = 7.2 MA. Current in units of 1.35 kA, frequency in units of 
Wo = 3 x 10 10 rad/s. 
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Figure 20. Periodic simulation AD, time history and frequency spectrum of the net cur- 
rent. Ii — 28.8 MA. Current in units of 1.35 kA, frequency in units of 
u) 0 — 3 x 10 10 rad/s. 
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Figure 21. Particle distribution phase-space plot of axial momentum (Pi = P„l) versus 
radial position (X2, cm). Left-hand plot is sheath electrons, right-hand plot is 
beam electrons, (a) and (b) 1.8 MA run AB at ct — 400 cm. (c) and (d) 7.2 MA 
run AC at ct = 500 cm. (e) and (f) 28.8 MA run AD at ct = 200 cm. 
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3.0 LAMINAR FLOW EQUILIBRIUM MODEL 



We propose a laminar flow model (similar to a model that has been used for magnetic 
Insulation calculations 2 ) as a first simple description of the electron equilibrium In the 
plasmold. By “laminar” we mean that v r = 0, Nevertheless, we assume that all electrons 
now In the plasmold of radius r* were born at zero energy on a zero-potential (<^ = 0) 
surface coinciding with a conducting pipe at radius > r&. Energy conservation thon 
gives the usual relativistic factor 

7 - 1 = c</>/md 2 (1) 

where e (> 0) and m are the electron charge and mass. The ion-beam density n,' and the 
axial velocity v M < are assumed to be free parameters Independent of radius for r < r> and 
v,i = 0 for r > r&. Conservation of axial canonical momentum for the plasmold electrons 
gives 



P M s 7 mu, - eA, — 0 



( 2 ) 



where we have taken the vector potential A, = 0 at the <t> = 0 surface. 



The radial electric field is defined In terms of the potential by 



E r = 



d(f> me 2 c<7 



dr 



e dr 



(3) 



and is given in terms of the plasma sources by Poisson’s equation 



Id/ d(f>\ _ e (n e - n<) 
r dr \ dr ) € 0 



( 4 ) 



Thus, using Eqs. 3 and 4 the electron density is 



eomc 2 1 d ( d^\ 
n ' = —rd-rVTr) +ni ' 



( 5 ) 
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The azimuthal magnetic held la defined in terms of the vector potential from Eq, 2 




m d 
e dr 



('m) . 



( 0 ) 



Combined with E r from Eq, 3 this gives the radial force-balance equation E r - v,Bt = 0, 
which gives E r w cBt In agreement with the ISIS simulation results. Ampere’s law relates 
Bt to the plasma sources by 



-■j- (rB f ) = n 0 e - n,v t ) . 
r ar 

Substituting Bt from Eq. 6 into Eq. 7 gives 



(?) 



1 d 
r dr 




m 



(n,f„ - n.f.) = 0 . 



( 8 ) 



Substituting the electron density n, from Eq. 5 into Eq, 8 gives a second-order radial 
nonlinear differential equation for the electron velocity v M \ 



\ d_ 
r dr 





(9) 



For simplicity we analyze the nonrelativistic limit 7 « 1, valid for v„ < c, which 
linearizes the equation. We also ignore cylindrical effects, valid for rwt/c » 1. Defining 
y = v M - Vg( and x = (r — r^)w(,/c, Eq. 9 reduces to 

TT? = y (x<0) (10) 



which has the general solution y = cjc* + cje - *. To keep y finite for large negative x 
(r < r/,), we require a — 0. Outside the beam, v M i — 0 and »,• may be different than Inside 
the beam. Defining X 2 = n,(r > ri)/n,(r < rj,), the electron velocity is 



d 2 v M 



= A 5 u r 



( x > 0) 



(H) 
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which haa the general solution v t =» Osfl A " + c<<T A *. We assume (r w - r^Xui^jo > 1 (many 
plasma skin depths between beam and wall), To keep v g finite for large positive x, we 
require 03 = 0 , Thus, the velocity Is described by 



v$ 



v„i + cie* (x < 0 ) 

i 

c<e" A " (x > 0 ) 



( 12 ) 



Because the potential (f> and Its derivative, the electric field, must be continuous at the 
beam edge x = 0 , this defines the two constants and gives 



v*_ 

Vgi 



1 “ ( I<o » r<r 0 
(x> 0 ,r>r 6 ) 



(13) 



Thus, at the x =■ 0 (r = rt) the velocity is 

(r « r t ) __ 1 

v,i 1 + A 



( 14 ) 



which Is small If X >■ 1 , while v, —*■ v,< for X — + 0 , 

Several relationships follow from this analysis. From Eq, 6 , the magnetic field at the 
beam edge Is 



„ , , mv t i X . , 

fl,( r = r ,) = __-_ u , ( 15 ) 

which vanishes for A -+ 0 and approaches a constant for A » 1 . The net current In the 
beam Is 



lntt (r =* ri,) 



~B,( r = n) = iiUfc 
Mo * 



X 

1 + A 



{Wb/c) 



(16) 



where I a 3 47rmc/e/,t 0 = 17 kA. Using the relation 




where /* Is the Ion beam current, Eq, 10 reduces to 





( 18 ) 



The azimuthal magnetic Held exerts a confining pressure fl|f(rj)/2Mo on the Ion beam 
which must be countered In equilibrium by a finite lon-beam pressure P{ = n<T<. Using 
Eq. 13 and equating these two pressures, the ion beam density cancels from the equation, 
leaving only an expression for the Ion temperature 

r < - I™'?'* (iTa)’ ' (19) 



This Implies that the Ion beam Is generally quite hot In equilibrium with 5=3 

A/(l + A). 

These theoretical predictions are specific and in some cases can be compared with 
existing ISIS simulations. The square of the net current / na « , divided by the beam current 
h, Is plotted versus It, for the simulation results In Fig. 17, and the result Is Indeed constant 
at approximately 5 kA Independent of beam current. For all of our plasmold simulations 
(here and in Ref. 1), we used an ion drift velocity of /?„< = 0.62. For all but the earliest 
simulations In Ref. 1, we use a plasma density between the wall and beam equal to five 
times the beam density, or A - V&, Therefore, Eq, 18 gives a predicted value of 5.03 kA 
for the constant of proportionality between and /*. Moreover, Eq. 19 gives a 

predicted ion-beam temperature for equilibrium of = 46.9 kV. In the ISIS simulations 
In Ref. 1, we varied TJ and empirically found that 43 kV gave a good equilibrium for beam 
currents of 22.5 kA, 112,5 kA, and 452 kA. In our present research reported here, this 
same value of 43 kV was used for the ion-beam temperature in all of the simulations and no 
noticeable contraction or expansion of the Ion beam was observed. These two agreements 
are excellent, particularly when the wide range of beam currents U = 22.5 kA-28,8 MA is 
considered. 




The radial profile of the electron axial velocity In Eq. 13 also can be compared with 
simulation, In all of the MRC ISIS simulations, A wan large (usually A 2 = 5), The 
description of the electron velocity v t In Eq, 13 Is consistent with the simulation results In 
that v„/v t i was small at the beam edge x =a 0 and Increased rapidly over a few beam skin 
depths e/uih until v M « In the beam Interior, The near equality of E r and oBf In the 
simulations also Id consistent with the laminar flow model used here, 

Several scaling predictions remain to be verified, One trend apparent from Eq, 13, 
but not yet verified by computer simulation, Is that V, at the beam edge increases as A 
decreases and v,(r =» n,) -+ as A — + 0, From Eq, 15, the magnetic field at the beam 
edge Is proportional to A for A < 1, and saturates at a constant for A > 1. The constant 
of proportionality /„ between ii, (r/>) and /j, 

/.3/ a A,(~J , (20) 

Is linearly proportional to the Ion beam axial velocity /?-< which was never varied from 0,62 
In our series of simulations. There la also a strong dependence of /„ on A, For A < 1, 
/„ Is proportional to A 2 — that Is, to the ratio of exte rior plas ma density to beam density. 
Combined with Eq. 18, this predicts (r^) r > r*), or to the number of exterior 

plasma skin depths across the beam. For A > 1, as In the simulations, (r*,) oc 
or to the number of beam skin depths across the beam. The ion-beam temperature for 
equilibrium la extremely sensitive to parameters, It Is proportional to and, for A < 1, 
to A 2 . For A > 1, It saturates at (T<) mU{ — 1/2 mvjh, These changes In scaling should be 
verified by simulation. 




4.0 CONCLUSIONS 



Increasing the beam current to 28.8 MA does not appear to strongly affect the equi- 
librium scaling discussed In our earlier reports. While there is some Indication that the 
current sheath thickness Is now not decreasing as rapidly at higher beam current as at 
lower values, the net current still scales as I nt t ~ AVom* The beam temperature required 
for equilibrium appears to remain constant at 43 kV independent of beam current. 

To try to shed light on the observed scalings In the ISIS code results, we completed a 
laminar flow model of the electrons. Specific predictions were obtained for radial profiles 
as well as for the net current and the lon-beam temperature necessary for equilibrium. In 
all cases, the predictions agreed very well with the ISIS simulations, and new parameter 
scalings were brought out which should be verified. 




5.0 FUTURE WORK 



There are many beam and background plasma parameters that affect the plasmold 
equilibrium. In spite of having carried out extensive and complex simulations over the last 
two years, some of these parameters have remained fixed. The laminar flow model pre- 
sented here Indicates several interesting parameter scalings that, if verified, would further 
our confidence in this model. We could continue simulations on to higher currents, but we 
believe that this scaling is adequately verified over the present range of h = 22.5 kA- 
28.8 MA. For instance, by moving the metal “core” boundary to larger radius (l.e., 
R = 32 cm) and limiting the plasma sheath thickness to four Debye lengths we could 
probably run a 115 MA case to t = 100 (3.3 ns) In about 15 CPU hours, but we ex- 
pect little new Information to result. A set of 7 MA cases allcwing variation in beam ion 
temperature and velocity, sheath electron temperature, and sheath thickness and density 
would indicate how sensitive the high-current equilibrium is to these parameters and would 
further test the laminar flow model in new areas. 

If an analytic model can be found that adequately predicts the wavelength of the 
fastest growing two-stream-like instability we should be able to resolve one wavelength in 
the axial direction sufficiently to determine a growth rate at moderately high (~ 2 MA) 
current. This would provide a useful data point In the regime where the “mixing” region 
thickness Is much less than the beam radius and whether or not it stays that way with the 
instability present. 
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APPENDIX A 

LOW-MODERATE BEAM-CURRENT SIMULATION RESULTS AT AFWL 



Following is an excerpt from the quarterly progress report for September 1988- April 

1989. 

A.i SIMULATION RESULTS 

To reduce computation time (by a factor of ~ 40) we have used a “nearly 1-D” 
model to generate equilibria. Cases for plasmoid propagation with ion currents ranging 
from 22 kA to 1.8 MA have been run thus far. This corresponds to the number of skin 
depths across the beam, wje^/c, in the range 2.9 to 26.0, respectively. During the course 
of simulating higher currents, we noted that noise levels in the diagnostics increased with 
the current to a point that made interpretation very difficult by 1.8 MA. Increasing 
the current further to > 7 MA would produce nearly unreadable results. We, therefore, 
decided to try reducing the noise to acceptable levels. A series of runs (on the AFWL 
CRAY) were done changing the number of grid points, the “current smoothing” and/or 
the number of particles-per-cell in the axial direction (XI, the dimension that is being 
approximately eliminated) and in the radial direction (X2). Most of these tests were done 
at 22 kA where the run time was short but the noise level still easily observable. The 
“current smoothing” namelist parameters for XI and X2 (NSMl, NSM2) and the number 
of particles-per-cell were adjusted separately or concurrently. We conclude that having 
much more than 10 particles per cell and extensive smoothing with NSMl - 1, NSM2 = 3 
reduced the noise considerably with only two cells (three grid points) in XI. This set of 
parameters increased the running time by about a factor of two to four over the original 
case, which still extrapolates to a run time for a 7 MA case of about 8 hours, a high but 
acceptable cost. We then set up a run at a current of 7.2 MA which had run on the AFWL 
CRAY to about one-fourth of the equilibration time when the “free” SDI computer charge 
account was shut down due to lack of funds. Results at that point were not close enough 
to equilibrium to provide meaningful values. (Our expectation of further funding of this 
account has not materialized after several months.) The following figures describe the 
results of our efforts on the other cases to date. 




A.l.i LOW CURRENT: h = 22.5 kA (e.g„ r h = 49 cm and n t = 10 9 cm" 3 ) 

We begin our Series of runs using the particle-in-cell code, ISIS 5.0, developed at Los 
Alamos National Laboratory. Results and parameters used in this series can be found in 
Table A-l. Since we were running ISIS at the Weapons Laboratory for the first time, we 
felt it was necessary to rerun a previous case to insure the validity of future results. We 
chose LANL run AM (7* = 22.5 kA) as our comparison case. 1 Changing the run ID to 
A A (Table A-l), we reran the simulation and got very good agreement with the previous 
LANL case AM (Figs. A-la-A-3a). Since the amount of smoothing used affects the codes’ 
run time, run AB was built by turning-off the smoothing in XI and X2 (NSMl = 0 and 
NSM2 = 0). The results from run AB were very noisy and considered unusable (Figs. A- 
lb-A-3b). Although the results from run AB were not good, it did show us that smoothing 
in the radial direction was crucial and had to be used to reduce noise. Run AL was created 
from run AA by reducing the number of axial cells to NXl = 2 and turning-off the axial 
smoothing (NSMl = 0), but setting the radial smoothing to NSM2 = 3. While the results 
were a bit noisy (Figs. A-lc-A-3c), the run time improved by a factor of 34 as compared 
to run AA. Some of the noise in run AL might be attributed to the fact that only one 
slice was plotted for the field plots, while five slices were plotted for runs AA and AB. In 
Fig. A-3, run AL had only two cells for the axial integration. Runs AA, AB and AL all 
ran to ct = 3000 cm and appeared to be in a steady state (see Fig. A-16). A time-average 
current profile versus radius for the I\, = 22.5 kA cases can be found in Fig. A-4. The 
three runs exhibit nearly identical time-averaged behavior. 

Because of the fast running time of run AL, it was felt a parameter study could be 
done quickly and efficiently using this case. Nine runs were created (AC-AK) by varying 
different parameters and running each case to ct = 400 cm (Table A-l). The results from 
these runs show little or no change versus run AL at this early time and were not included 
in this report. 

A. 1.2 MODERATE-CURRENT: = 452 kA (e.g, = 49 cm and rtj, = 2 x IQ 10 cm' 3 ) 

For the 7j = 452 kA runs in this series, we used run AP (previously ran at LANL) as 
our base case, 1 All runs in this series (BA1, BC-BG, BI) used two axial cells, except for run 
BA2, which used four axial cells. Simulation runs in this series were carried out to different 



times. While most runs were carried out to either ct = 1600 cm or ct = 2000 cm, several 
runs ran only to ct = 400 cm. Figures A-5a-A-7h show B 0 versus r, E r versus r, and 
axially-integrated J M versus r for comparing runs (BA1-BA2, BC-BG, BI) at ct = 400 cm, 
while Figs. A-8a-A-10f show the same type of plots comparing runs (BA1-BA2, BC-BF) 
at later times (ct = 1600 cm and ct = 2000 cm). Results and parameters for these runs 
can be found in Table A-l. For this current, a steady state for the net current is reached 
at about ct « 800 cm (see Fig. A-17). 

Run BAl was the first run in the series and was used to compare with LANL run 
AP. 1 The results of BAl were much noisier than AP (Figs. A-5a-A-10a). Run BA2 was 
built from BAl by changing IQ = 2 to IQ = 1 (the IQ parameter is the integer iteration 
count for subcycling the particles; higher numbers resolve cyclotron motion better) and 
increasing the number of axial cells to NXl = 4. Also the radial smoothing was turned- 
off (NSM2 = 0). By comparing BAl (Figs. A-5a-A-10a) to BA2 (Figs. A-5b-A-10b) it 
shows the importance of having IQ > 2 and the need for radial smoothing (NSM2 > 0). 
Run BC was setup by increasing the number of axial particles per cell. The results from 
run BC (Figs. A-5c-A-10c) are quieter than run BAl (Figs. A-5a-A-10a) and shows 
that increasing the number of axial particles per cell will help reduce noise. Run BD 
was created by modifying run BC to become a “hollow-beam” case, that is no particles 
were allowed in the first row of cells along the axis. The results of run BD (Figs. A-5d- 
A-lOd) show little or no improvement using the hollow-beam. Run BE was generated 
by reducing the axial and radial cell size from 0.2 cm to 0.15 cm. The results from run 
BE (Figs. A-5e-A-10e) are quieter than run BC (Figs. A-5c-A-10c) and indicate the 
need for a reduced cell size. Most runs to this point used radial smoothing (NSM2 = 3), 
but had the axial smoothing turned-off (NSM1 = 0). Run BF (Figs. A-5f-A-10f) was 
built from run BE (Figs. A-5e-A-10e) by setting NSM1 = 1, and the results showed a 
slight noiSe reduction by using axial smoothing. Feeling that the best way to continue to 
reduce noise was to increase the number of axial and radial particles per cell, run BG was 
created in this manner. By comparing the results of run BG (Figs. A-5g-A-7g) to run BC 
(Figs. A-5c-A-7c), it is apparent that increasing the number of particles per cell does 
reduce noise. Run BI (Figs. A.-5h-A-7h) was built from run BG (Figs. A-5g~A-7g) by 
decreasing the number of axial particles per cell by one half. The results showed an increase 
in noise and further proves the importance of having a sufficiently large number of axial 
particles per cell. Comparing run BI (Figs. A--5h-A-7h) with run BAl (Figs. A-5a through 




A-7a) shows the importance of having a sufficiently large number of radial particles per 
cell. Figure A-ll showB time-average current profiles for runs (BA1-BA2, BC-BF) which 
ran to times > 1600. Again, all the runs exhibit nearly identical time-averaged behavior, 
except for run BA2. 

A.l. 3 HIGH CURRENT; h = 1808 kA (e,g„ n = 49 cm and n t = 8x IQ 10 cm' 3 ) 

Only one run was made at /* = 1808 kA before the Weapons Laboratory account 
expired. Run BB (Figs. A-12a-A-12f) ran out to ct — 2000 cm. The net current appeared 
to reach a steady state at about ct « 400 cm (see Fig. A-18). This run had IQ = 1 and 
too few particles per cells radially, and the results looked similar to runs BAl (Figs. A-5a- 
A-lOa) and BA2 (Figs. A-5b-A-10b) in quality. Because of the extreme noise, the only 
results that were of much use was the time average current profile (Fig. A-13) taken at 
ct = 2000 cm. 

A. 1.4 CONCLUSION 



It is apparent that several parameters have a great influence on reducing the amount 
of noise in a run. Having the axial and radial smoothing (NSMl and NSM2) turned-on is 
crucial to noise reduction, and also IQ (the integer iteration count for subcycling particles) 
must be set to a value > 2. The reduction in cell size was important in improving the 
simulation results. Having a sufficiently large number of particles per cell appears to have 
had the most noticeable effect of all. 

The scaling of the current radial-profile and the net current is summarized in two 
figures. Figure A-14 shows the time average current profile for different values of the beam 
current, and Fig. A-15 shows the net current as it scales to the square root of the beam 
current. 

Figures A-16a-A-16b show a good comparison of the effect that the number of 
axial cells has on the noise in a simulation. Figure A-16a is a time history and frequency 
spectrum of the net current for run AA with 129 axial cells with a beam current of 22.5 kA. 
Figure A-16b is a time history and frequency spectrum of the net current for run AL with 
two axial cells with a beam current of 22.5 kA. After examination one can see that the 




number of axial cells has only a slight effect on the noise. Figures A-17 and A-18 show 
time history and frequency spectrum of the net current for beam currents of 452 kA (run 
BF) and 1808 kA (run BB), These figures show that the noise increases as the beam 
current increases unless more cells and particles are used. The net current values used 
In Figs. A-14-A-15 were taken from the time history plots shown In Figs. A-16-A-18. 
Figures A-16-A--18 also show that the net current reaches a steady state in a time of order 
ct/rt, - 200 c/ub,ri,. Thus, the higher current runs may not have to be run as long. As 
a rough estimate, the number of CRAY hours to reach steady state appears to scale as 
^rtm (hr) '•sV 2(wi,,r(,/l3c) (NPC/40), where NPC is the number of particles per cell. 

By taking an overall view of the runs done in h — 452 kA series, run BG (Figs. A- 
5g-A-7g) was giving the best results followed by run BI (Figs. A-5h-A-7h) at early times 
( ct — 400 cm). At later times ( ct > 1600 cm) run BF (Figs. A--5f-A-10f) gave very 
good results. Again the major factor in all cases that returned reasonable results was that: 
(a) smoothing was turned-on, (b) IQ was set to at least two, (c) a sufficiently large number 
of particles per cell was used (> 40). 

In order to achieve results comparable with run BG (Figs. A-5g-A-7g) for 1808 kA, 
run BB (Figs. A-12a-A-12f) would have to have IQ = 2, DXl = DX2 = 0.1 cm, six 
particles per cell axially and seven radially (42 total). A simulation with these parameter- 
settings would take about 4 hours of Cray time to execute out to a time of ct = 400 cm. 




TABLE A-l. SIMULATION RESULTS AND PARAMETERS. 







• l 



6 



# 



# 




CURRENT (kA>) 



LLNL PLASMOID PROPAGATION 




10 20 30 40 50 60 

RADUS (cm) 



RUM AA 

1 — 

RUN AB 



RUM AL 



Figure A-4. Time averaged current profile in steady state for runs at h = 22.5 kA. 
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Figure A-15. Net current as it scales to the square root of beam current. 








Figure A-16, Periodic simulation results comparing time history and frequency spec- 
trum of the net current for two runs with a different number of axial 
cells, h = 22.5 kA. Current in units of 1.35 kA, frequency in units of 
Wo = 3 x 10 10 rad/s. (a) run A A with 129 axial cells, (b) run AL with two 
axial cells. 
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Figure A-17. Periodic simulation BF, time history and frequency spectrum of the net 
current. = 452 kA. Current In units of 1.35 kA, frequency In units of 
w 0 = 3 x 10 10 rad/s, 
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Figure A-18. Periodic simulation BB, time history and frequency spectrum of the net 
current, h = 1808 kA. Current In units of 1,35 kA, frequency In units of 
w 0 = 3 X 10 10 rad/s. 








